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ABSTRACT 

Aims. We report on a detailed study of the optical afterglow of GRB 061 121 with our original time-series photometric data. In conjunction with 
X-ray observations, we discuss the origin of its optical and X-ray afterglows. 

Methods. We observed the optical afterglow of Swift burst GRB 061121 with the Kanata 1.5-m telescope at Higashi-Hiroshima Observatory. 
Our observation covers a period just after an X-ray plateau phase. We also performed deep imaging with the Subaru telescope in 2010 in order to 
estimate the contamination of the host galaxy. 

Results. In the light curve, we find that the optical afterglow also exhibited a break as in the X-ray afterglow. However, our observation suggests a 
possible hump structure or a flattening period before the optical break in the light curve. There is no sign of such a hump in the X-ray light curve. 
Conclusions. This implies that the emitting region of optical was distinct from that of X-rays. The hump in the optical light curve was possibly 
caused by the passage of the typical frequency of synchrotron emission from another forward shock distinct from the early afterglow. The observed 
decay and spectral indices are inconsistent with the standard synchrotron-shock model. Hence, the observation requires a change in microphysical 
parameters in the shock region or a prior activity of the central engine. Alternatively, the emission during the shallow decay phase may be a 
composition of two forward shock emissions, as indicated by the hump structure in the light curve. 
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1. Introduction 

Gamma-ray bursts (GRBs) and their afterglows are widely be- 
lieved to be emission from relativistically ex panding shells (e.g., 
Zhang & Meszaros 20041 iMeszaros 1 12006). GRBs, or prompt 
emissions are considered to arise from internal shocks caused 
by collisions between the shells. After the collisions, the shell 
keeps expanding and generates an external shock colliding with 
the interstellar medium. As a result, synchrotron emission from 
the shocked region is observed as afterglows. This synchrotron- 
shock model successfully reproduces the observed temporal evo- 
lution of spectral energy distributions (SEDs) of late afterglows. 
According to this model, the flux of the synchrotron emission 
from afterglows is described with a power-law form, that is, 
f v (t) oc t~ a v ~P, where a and /? are a de cay index and a spec- 
tral slope, respectively (Sar i et all Il998l) . A jet geometry was 
suggested by an achromatic break observed in light curves of 
afterglows (e.g., Rhoads, 1997; Sari et al., 1999). This has re- 
cently been called into question because chromatic breaks were 
detec ted at the time when th e jet scenario predicts achromatic 
ones dWillingale et al.ll2007l) . 

Owing to quick identifications and notifications of GRBs by 
the Swift satellite, the number of observ ations of early afterg lows 
has been increasing in all wavelengths (Geh rels et all 120041) . X- 
ray light curves of early afterglows, in particular, turned out to 
have more complicate profiles than those previously expected 
from the standard synchrotron-shock model. Although a simple 
power-law decay was expected in the standard model, the early 
X-ray light curves actually consist of three stages with different 
decay indices; the initial steep decay (a ~ 3 - 5), the shallow 
decay (a ~ 0.5 - 1.0), and the normal decay phases (a ~ 1) 
dNousek et alJ2 006; O' Brien et ai1l2006l) . While the steep decay 



phase is likely a high-latit u de emission of the prompt emission 
( Kumar. Panaitescul 120001 ; lYamazaki et alj 1200 6; Liang et al. 
2006; Zhan g et al. I l2006h . the origin of the shallow decay phase 
is currently unknown. Several models have been proposed for 
this phase, for example, the late energy injection into the 
shocked region, the time-d ependent microphysic s in the shock, 
or prior outflow emission (Sari, Meszaros 2000; Nousek et al 



2006; Zhang, Meszaros 2001; Panaitescu et al. 2006; Toma et al 
l20061:lloka et al.ll2006l:lDado et alJl200dl Yamazaki ft2009l) . 

Each of those models predicts distinct behaviors of SED 
variations in early afterglows. Simultaneous multiwavelengths 
observations are required to provi de crucial clues on the na- 
ture of the early afterglow phase. Panaitescu et al. (2006|) re- 
ported that the early break after the shallow decay pha se is chro- 
matic on the basis of 6 afterglows. lYost et all J2007) addition- 
ally reporte d that GRB 05 11 09 A a lso exhibited a clear chro- 
matic break. Panait escu et all (120061) proposed that the observed 
light curves require the temporal evolution of microphysical pa- 
rameters in the emitting region of early and late afterglows. On 
the other hand, some afterglows apparently exhibited achromatic 
breaks after the X-ray shallow decay phase. Krii hler et alj ([2009) 
reported on optical — IR and X-ray light cur ves of GRB 080710 , 
in which an achromatic break was observed. Blus tin et al.l ((2006) 
reported another example of a possible achromatic break in 
GRB 050525 A, while X-ray flares make difficult to accurately 
determine a break time (also see. iKlotz et al.l 12005). In some 
of past cases, optical observations were too sparse to determine 
break times and to catch the detailed behavior on either side of 
the breaks. We definitely need new observations in which break 
times can be determined accurately both in X-ray and optical 
light curves. 
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GRB 061121was detected by the Swift Burst Alert Telescop e 
(BAT) at 15:22:29 (UT) 21 November 2006 (iPage et al.l 12006). 
Swift also reported the discovery of a bright optical afterglow 
with Ultraviolet/Optical Telescope (UVOT), which was soon 
confir med at 14.9 magQby ground-based telescopes dYost et all 
2006). Its redshift was estimated to be z = 1.314 by spectro - 
scopic observations of the optical afterglow dBloom et al.l l2006). 
This bright burst is a typical GRB followin g a well-known em - 
pirical relationship between E p and E lso dAmati et al.l [2006). 
The early X-ray light curve has several breaks as other systems 
observed in the Swift era dNousek et al.l 120061) . Prompt onsets 
of multi wavelengths observations for GRB 061121 provided a 



Table 1. Results of our photometric observation. 



unique opportunity to study the tempor al evolution of X-ray 

:al.l l2006t 



and optical afterglows (lYost et al.l I20p6t iMelandri eti . 
lUemura et al. 20061). IPage et all d2007l) reported on multiwave- 
length data during the prompt and afterglow phase of this GRB. 
According to them, both X-ray and optical flux monotonically 
decayed, which can be described with an early exponential rise 
followed by a power-law decay phase. 

Here we report on our optical and infrared observations using 
the Kanata 1.5-m telescope. Our continuous time-series obser- 
vations enabled us to reveal the optical behavior near the X-ray 
shallow decay phase. We describe the details of our observations 
in section 2. Combined with other published data, we report the 
temporal evolution of the optical and X-ray afterglows in sec- 
tion 3. In section 4, we discuss the nature of the variations in 
the light curves using the synchrotron-shock model. Finally, we 
summarize our results in section 5. 



2. Observation and Data Analysis 

2.1. Optical observations 

Our observation started at 16:37 (UT) 21 November 2006, 
4.6 x 10 3 s after the GRB trigger time, and ended at 19:57 (UT). 
The observation was performed with TRISPEC attached to the 
Kanata 1.5-m telescope at Higashi-Hiroshima Observatory of 
Hiroshima University. TRISPEC is a simultaneous imager and 
spectrograph with po larimetry covering bo th optical and near- 
infrared wavelengths (Watan abe et al.l l2005). We used the imag- 
ing mode for the observation of GRB 061121 and obtained 77 
sets of R c , J, and K s band images. The exposure time of a R c - 
band image was 123 s. During the 123 s exposure, short expo- 
sures of a few seconds were taken for NIR arrays, and yielded 
net exposures of 120 and 96 s for each J- and A^-band image, 
respectively. 

The central wavelength of the TRISPEC's R c system is ~ 
620 nm, slightly shifted from the standard one (= 645 nm). The 
difference in magnitude between these systems is expected to be 
0.008 mag when a power law spectrum with a spectral index of 
1 .0 is assumed. In the following discussion, we neglect this small 
difference. 

We show an example of the obtained images in the right 
panel of figure Q] We also show the same field in the Second 
Palomer Sky Survey (POSS2) in the left panel for compari- 
son. The afterglow is the object marked with the black bars. 
After making dark-subtracted and flat-fielded images, we ob- 
tained magnitudes of the afterglow and comparison stars us- 
ing a Java-based PSF photometry package. For a comparison 
star, we used USNO-B1.0 0768-0239968 (R.A.= 9 /! 48 m 54 s .78, 
Dec.= -13°1'17".9; R c = 18.02). The comparison star was 
constant within 0.02 mag during our observation, checked 



Time (s)* 


R" mag. 


R c mag 1 . 


error 


N+ 


4699 


18.16 


18.17 


0.03 


3 


5283 


18.23 


18.24 


0.15 


4 


ozol 


1 n in 

15.4U 


1 o /to 
15.42 


n no 
O.OO 


I 


6675 


18.53 


18.55 


0.09 


5 


7501 


18.63 


18.66 


0.09 


6 


8427 


18.74 


18.76 


0.18 


7 


9466 


18.77 


18.80 


0.10 


8 


10632 


18.92 


18.95 


0.10 


9 


11941 


19.18 


19.22 


0.09 


10 


13409 


19.49 


19.55 


0.23 


11 


15056 


19.10 


19.14 


0.21 


12 



* Time since the GRB trigger. 
** Raw R c magnitude, 
t R c magnitude (host corrected). 
t Number of images in each bin 



by USNO-B1.0 0767-0229365 (R.A .= 09 ,, 49"W.080, Dec.= 
-13°13'22".21). Using neighbor USNO stars, we checked sys- 
tematic errors of magnitudes depending on comparison stars, 
and found that it is smaller than 0.2 mag. The c ompari - 



Mag show Vega magnitude in this paper. 



son star is the same a s that u sed in Haloern et al.l (120 06a). 
lHalpern and Armstrong! (l2006al) . and lHalp ern and Armstrong 
(2006b) which present observations in a late stage of the after- 
glow. In the following section, we performed an analysis of our 
light curve in conjunction with those late-time observations. 

Additionally, tableQ]contains R c -magnitudes obtained by our 
optical observations. In this table, the magnitudes are averages in 
equally spaced bins in the logarithmic scale of the time. While 
we obtained J and K s band images using TRISPEC, IR after- 
glows were not significantly detected. Typical 3-sigma upper 
limits of each frame are 16.0 and 13.7 mag in J and K s bands, 
respectively. 

For deep photometry of the host galaxy component, we 
obtained /? c -band images with the 8.2-m Subaru Telescope 
and the Faint Object Camera and Spectrograph (FOCAS; 
iKashikawa etai1l2002l) on 2010 May 7 (UT). The total expo- 
sure time was 240 s. We can easily recognize the host galaxy 
as a point source at the GRB afterglow position in the obtained 
image. Using the same comparison star as mentioned above, we 
derived the magnitude of the host galaxy to be R c - 22.99+0.03. 

2.2. Data analysis of the Swift data 

We analyzed public data of GRB 061121 observed with X- 
Ray Telescope (XRT) and UVOT on Swift. We processed the 
XRT all orbits of data, adopting the standard screening with 
the XRT pipeline FTOOL xrtpipeline (Version: 0.10.3). We ex- 
tracted light curve and spectra with a rectangular 40x20-pixel 
region for the Windowed Timing (WT) mode, and 40-pixel ra- 
dius region for the Photon Counting (PC) mode from the source 
position, respectively. The background was also extracted from 
40x20-pixel region for the WT mode, and 40-pixel radius region 
for the PC mode, far from the source, respectively. While begin- 
ning of GRB for WT mode data and PC mode data, we found 
that the count rate is high enough to cause the pile-up effect, 
and we adopted the s tand ard pile-up correction as described by 
iRomano et all d2006l) and lVaughan et al.l d2006l) . 
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Fig. 1. Optical images of the field of GRB 061121 in POSS2 (left panel) and observed with the Kanata 1.5-m telescope (right 
panel). The field of view is 4' x 3' and the top is north. The afterglow is an object marked with the black bars in the right panel. 



In the following section, the unit of time is set to be seconds 
from the GRB trigger. Optical and X-ray parameters are indi- 
cated by subscripts of "O" and "X", respectively. 



3. Results 

3. 1 . Optical and X-ray light curves 

We show the X-ray and optical light curves of GRB 061121 
in figure [2] Our optical observations and X-ray obser- 
vations by XRT are indicated by the filled circles and 
crosses, respectively. The open circles and squares are 
optical observations with R c , V/"White light" (UVOT ) 
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2006; 



bands reported to GCN, respectivel y (Page et al 
Yost et all [20061: iMelandri et all 12006 1 : ISonoda et al. 
Marshall et al J 120061: iGolenetskii et al J 120061: iHafoern et al] 
2006 bt ICenkoll2006t Habern et alj|2006al: lEfimov et al.l l2006a; 



Efim ov et al] l2006bl) . About 



Halpern and ArmstrongT 2*006allbt 
the optical flux, the contribution of the host galaxy is subtracted. 
The flux density of the host galaxy is 2. 14 fiJy in R c , which was 
estimated based on our Subaru observation as described in sec- 
tion 2. The flux dens ity was corrected for the Galactic extinction 
of E(B - V) = 0.04 dSchlegel et al.lll998l) . The absolute magni- 
tude of the host g alaxy is -21 . 86 m ag in R c . 

According to Page et aD d2007l) . the X-ray light curve of 
GRB 061 121 is divided into 4 phases depending on their decay 
indices. In this paper, we follow their definition of the phases 
for the X-ray light curves, that is, an initial flare, a plateau 
(axi = 0.38 + 0.06), a shallow decay (a x2 = 1.07 ± 0.05), and 
a normal decay phases (ffx3 = 1.53 + 0.03). The errors of these 
parameters, as well as other parameters given in this paper, rep- 
resent 1-cr. 

In the optical light curve, we can see a possible flare in a 
very early phase at t — 76 s. This implies that the optical flux 
may be associate d with the promp t emission in X-rays and y- 
rays in this phase (Page et al.Ll2007T) . In the X-ray plateau phase, 
the optical light curve can be described with a simple power-law 
decay. Using V-band observations by UVOT from 240 to 2000 s, 
we calculated the power-law decay index to be ao\ = 0.72±0.08. 

During the subsequent shallow decay phase, our observation 
revealed a monotonic fading of the optical afterglow. The light 
curve can be described with a simple power- law having a decay 
index of a 02 = 0.96 + 0.06. lPage et alJ (120071) reported a V-band 
decay index of 0.66 + 0.04 from the onset of the fading to a 
break at ~ 2.5 x 10 4 s. This decay index was estimated based on 



the exponential-to-power-law model, which assumes a mono- 
tonic fading during the fading stage of the afterglow. The ao2 
estimated from our time-series ph otometry is, how ever, signifi- 
cantly larger than that reported in Pag e et al.l d2007l) . We tried to 
fit a simple power-law model to the optical light curve from 240 s 
to 16 ks including our data. For the fitting, our V-band data was 
shifted to corresponding /? c -band magnitude. The V - R c of the 
afterglow was estimated from two almost simultaneous V- and 
/? c -band observations, that is, t ~ 280 and ~ 6300 s. The aver- 
age color of the afterglow is calculated to be V - R c = 0.05. The 
best-fitted parameters yield a chi-square/d.o.f of 76.3/15. This 
value is too high to conclude that the afterglow decayed with a 
simple power-law form from 240 s to 16 ks, and rather suggests 
that there is a sub-structure around the termination of the X-ray 
plateau phase. 

The optical light curve, then, exhibit another break around 
3 x 10 4 s, which is followed by the normal decay phase described 
with ao3 = 1 .58+0.03. By fitting a broken power-law model, we 
calculated a break time of 4. 6^3 x 10 4 s. About the last break 
from the shallow decay to the normal decay phase, no signifi- 
cant time lag is detected between the X-ray and optical breaks, 
while the errors of break times are quite large; 3.2^ x 10 4 s and 
4.6^ 4 3 x 10 4 s for the X-ray and optical break times, respectively. 
It is noteworthy that the optical decay index is almost same as 
the X-ray one after this break. According to the standard syn- 
chrotron shock model, this strongly indicates the passage of the 
cooling frequency of the synchrotron emission within the optical 
band at the break time. 

We searched possible correlations between X-ray and opti- 
cal short-term variations. We calculated cross-correlations using 
3 segments in which simultaneous optical and X-ray data are 
available. The optical data was divided into the following parts; 
i)4.8xl0 3 s < t < 7.5xl0 3 s,ii) l.OxlO 4 s < t < 1.4xl0 4 s, and 
iii) 1 .6 x 10 4 s < t < 1 .9 x 10 4 s. The resultant cross-correlations 
are shown in figure[3] The correlation functions are flat and show 
no prominent feature. We cannot detect any significant correla- 
tions between optical and X-ray short-term variations. 



3.2. Spectral energy distribution 

Figure|4]shows infrared-X-ray SEDs. The figure contains 5 pan- 
els in which simultaneous optical and X-ray observations are 
shown for the 5 phases. We fitted a power-law model with a sin- 
gle absorption component for the X-ray spectra. All 5 spectra 
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Fig. 2. Optical and X-ray light curves of GRB 061 121 afterglows. The abscissa denote the time since the GRB trigger in seconds. In 
the top panel, the ordinate denote the flux density in mJy for optical data and the count rate for X-ray data observed with XRT. Our 
observations and X-ray observations by XRT are indicated by the filled circles and crosses, respectively. Open circles and squares 
are optical observations with R c and V/White light (UVOT) bands reported to GCN or taken by UVOT, respectively. The solid lines 
are 95 % confidence regions of the best fitted power-law models for optical light curves. The labels, (a),(b),(c),(d), and (e), represent 
the time intervals for the SED analysis (see section 3.2). In the bottom panel, the ordinate denote the ratio of the observed flux 
density to the best-fitted model of the plateau phase. V band points are shifted by -0.05 to match the R c band points. 




1000 



100 



Optical lag ( 



Fig. 3. X-ray/optical cross correlations for short-term variations. The solid, dashed, and dotted lines were calculated from the data 
of 4.8 x 10 3 s < t < 7.5 x 10 3 s, 1.0 x 10 4 s < t < 1.4 x 10 4 s, and 1.6 x 10 4 s < t < 1.9 x 10 4 s, respectively. The upper and lower 
panels show correlation functions in long and short time scales, respectively. No significant correlation can be seen. 



in figure |4] can be described with absorption models with a hy- 
drogen column density of Nh = 2.2+ 0.15 xlO 21 crrr 2 in the 
observer's frame, which corresponds to Nh = 9.2 X 10 21 cm -2 in 
the rest-frame. The solid lines in the figure indicate the best fit- 
ted unabsorbed power-law component of X-ray spectra. As can 
be seen from the figure, spectral slopes /3 (f v oc v~@) were slightly 
larger in the plateau phase (panel b) and a phase just after the op- 



tical break (panel c) than those in later phases. Table [2] contains 
results of our best fitted parameters in each period. 

In the figure, the optical flux was corrected for the galac- 
tic and extragalactic extinctions. The correction for the extra- 
galactic extinction was performed with the relationship between 
the visual extinction A y and the hydr o gen c olumn density Nh 
for the "Q2" model in iMaiolino et al.1 ( feoOll) (N H /A V = 3.3 X 
10 21 cirT 2 ). We estimated the Nh from the best-fitted model of 
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Table 2. Best fit parameters for X-ray spectra 
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X-ray spectra, and set to Nh = 9.2 x 10 21 crrT 2 as mentioned 
above. The conversion from Ay to those in other bands was per- 
formed following the equations in lCardelli et alj d 19891) . 

For the correction of the extragalactic extinction, we used 
the "Q2" model because it provides the most plausible optical — 
IR SEDs, as shown below. According to the synchrotron-shock 
model, the spectral slope at the optical region should be Po = 
fix - 0.5 in the case of v m < vo < v c < vx, where v m and v c 
are typical and coolin g frequencies of the synchrotron emission 
from a forward shock ( Sarietalfl 998). In figure |4l we show the 
expected spectral slope of fio = fix - 0.5 with v c = 10 18 Hz, 
indicated by the dotted lines. The optical flux is required to 
be over this dotted lines in order to satisfy the condition ex- 
pected by the synchrotron-shock model. In addition to the "Q2" 
model, the figure also contains optical-IR points corrected with 
the Milky Way model (N H /A V = 1.6 X 10 21 crrT 2 ; open tri- 
angles) and w ith the "Ql" model (7Y H /A V = 6.7 x 10 21 cirT 2 ; 
open squares)) Mai olino et al.ll200lh . As can be seen in panel (c), 
a high Ay provided by the Milky Way model yields an unnat- 
urally sharp break between R c and /-bands. Corrected by the 
"Ql" model, the optical flux is too low to be interpreted by 
the synchrotron-shock model with fio = fix - 0.5. The mod- 
els for SMC (N H /A V = 1.5 x 10 22 cirT 2 ) and LMC (N H /A V = 
7.6 x 10 21 crrT 2 ) also yield further lower optical flux. Thus, the 
"Q2" model provides the best correction among those models. 

Near the peak of the prompt emission, as can be seen in 
panel (a) of figure |H the optical flux is much above the power- 
law component of X-rays. This indicates that the emission mech- 
anism or source of the optical emission is distinct from those of 
the prompt X-ray and y-ray emission. 

In panel (c) of figure [4] there is a difference between the 
optical flux and the spectrum extrapolated from the X-ray data 
(the solid line in the figure). The SED, hence, requires a spectral 
break between the optical and X-ray bands. This is consistent 
with the situation for the case of v m < vo < v c < vx- In the 
standard model, v c evolves with time, decreasing in case of con- 
stant density medium, increasing in case of wind medium. Since 
our findings privilege a constant density medium as discussed in 
subsection 4.2, we expect at one point that v c , decreasing with 
time, will cross the optical band. At that time, the optical and X- 
ray decay will become identical, and the SED will be compatible 
with a simple power law. This is the case in panel (e), suggesting 
that v c crossed the optical band around ( = 5x 10 4 s. 

4. Discussion 

4.1. A possible hump structure in the optical light curve 

Here, we discuss the substructure in the optical light curve 
around 5 ks. The behavior of the optical and X-ray light curves 
is unclear in the transition phase from the plateau to the shallow 
decay phases (1.5 x 10 3 s < t < 4.6 x 10 3 s). In figure|2] the 



solid lines denote 95 % confidence regions of the power-law de- 
cay model for each phase of the optical afterglow. As indicated 
by these lines, our observation suggests that the plateau phase is 
terminated by a flattening phase or a hump in the optical after- 
glow. 

We checked the significance of this hump with a correction 
of colors, because the best-fitted models were calculated with 
data taken with different bands; the V- and 7? c -band observations 
in the plateau and shallow decay phases, respectively. The lower 
panel of figure [2] shows the ratio of the observations to the best- 
fitted model of the plateau phase. We converted V-band obser- 
vations to our /? c -band ones, by adding V - R c — 0.05 to the 
V-magnitudes. The hump appears over a 2.2-sigma level even 
with the color correction. 

By contrast, there is no sign of such a hump in the X-ray 
light curve between those two phases. In addition, there is a large 
difference in the optical and X-ray decay indices in the plateau 
phase, compared with those in the subsequent shallow and nor- 
mal decay phases. The X-ray and optical light curves apparently 
exhibit different behaviors during the X-ray plateau phase and 
the transition phase to the shallow decay phase. 

4.2. Implication to the synchrotron-shock model 

In this subsection, we discuss the behavior of the afterglow of 
GRB 061 121 based o n the synchrotron- shock model. In the case 
ofv m < vo < v c < vv. lUrata et alJ d2007l) proposed a relation be- 
tween the decay indices of the X-ray and optical bands described 
as ax - co = 1 /4. In the plateau and shallow decay phases of 
GRB 061121, the a x - a is 0.34+ 0.10 and 0.11+ 0.08, re- 
spectively. The classical synchrotron-shock model, hence, fails 
to reproduce the observed light curves of GRB 061121 in the 
shallow decay phase. It can only marginally reproduce the light 
curv e in the plateau p hase. 

iPanaitescu et alj d2006l) generalized the formulae of the 
synchrotron-shock model by including the variations of the en- 
ergy (E) in the blast wave, the energy ratio for electrons and the 
magnetic field (e, and eg), and the ambient medium (re) in the fol- 
lowing form; E(> T) oc T e , eg k F~ 6 , e,- oc r _! , and n{r) oc r~ s . 
The decay indices of optical and X-ra y light curves ar e calcu - 
lated as in equations (9) and (10) in Panaitesc u et alj J2006). 
Using those formulae, we evaluate the presence of the energy 
injection (<? > 0) or the time variations of microphysical param- 
eters (b * or i + 0) for GRB 061121. 

We first assume 5 = 0, namely the uniform distribution of the 
interstellar medium. In the following examination, we calculate 
p from /J in each time-interval shown in table [2] For the plateau 
phase, assuming b — and i = 0, we find that the optical and X- 
ray light curves yield inconsistent e, that is, e — 1 .32 + 0.24 and 
4.05+ 0.37 calculated from aoi and axi, respectively. This in- 
consistency can be reconciled only when s takes a narrow range 
of s = 1 .22 + 0.01 and e takes an unnaturally large value (<? ~ 7). 
In the other case of s > 0, the inconsistency in e becomes more 
extreme. These results indicate that the observed light curves 
during the plateau phase cannot be reproduced only with the en- 
ergy injection scenario. Temporal variations of the microphysi- 
cal parameters are, hence, required. Assuming e — and 5 = 0, 
we obtain b = -2.10 + 0.36 and i = 2.04+ 0.12. The positive i 
implies a low efficiency of the energy for accelerating electrons 
in highly relativistic shocks. 

During the shallow decay phase after the optical hump, the 
decay indices changed to ao2 = 0.96 and axi = 1-07. Even in 
this phase, the condition of b — and i - yields different e 
values calculated from C02 (yielding e = 0.64+ 0.20) and axi 
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v (rest-frame; Hz) 



Fig. 4. Spectral energy distributions derived for the time intervals of (a) 6.7 x 10 1 s < t < 8.7 x 10 1 s, (b) 4.3 x 10 2 s < t < 5.4x 10 2 s, 
(c) 5.0 x 10 3 s < t < 6.0 x 10 3 s, (d) 8.0 x 10 3 s < t < 1.2 x 10 4 s, and (e) 6.5 x 10 4 s < t < 1.0 x 10 5 s. The time intervals of the 
SEDs for each panel are indicated in figure 2. The abscissa and ordinate denote the rest-frame frequency in Hz and the flux density 
in mJy, respectively. The filled circles are X-ray observations by XRT, optical observations reported to GCN, and our optical — IR 
observations. The IR observations just provide upper l imits of the flux. For the extragalactic extinction, the optical and IR data were 
corrected with the "Q2" model dMaiolino et alJl200Th . The open triangles and squares are data in which their c orrections for the 
extinctions were performed with the Milky Way model and the "Ql" model, respectively dGuidorzi et alJl2007h . We omit arrows 
indicating upper limits of the IR points of the open triangles and squares. The solid lines represent the best fitted power-law model 
for the X-ray spectra. The dotted lines are expected spectral slopes from the synchrotron-shock model. 



(e = 1.19+ 0.25). Assuming e = 0, we obtain b = -0.50+ 0.27 
and i = 0.75+ 0.12 for this phase. It is interesting to note that 
the absolute values of both b and i decreases from the plateau 
to shallow decay phases. From the early to late stages, the con- 
dition of the blast wave may resemble the classical picture in 
which no temporal variation in the mic rophysical parameters is 
required. Thus, the model proposed by Panaite scu et alj (2006) 
can explain the light curves in both the plateau and shallow de- 
cay phases by changing the microphysical parameters. 

An alternative model was proposed by lloka et aP d2006h . 
which consider prior activities before the main prompt emission. 
According to their model, the shallow decay of X-ray afterglows 
appears because a blast wave obtains additional energy by col- 
liding with prior ejecta without significant decelerating of shells. 
This model is possibly preferable for GRB 061 121, since it has 
a precursory 75 s before the peak of the main prompt emission 
(Page et alj|2007l) . The precursor may be a sign of the existence 
of the prior activity. 

lloka et alJ J2006h define the the prior ejected mass as a 
power-law form of y, that is, M(< y) oc y" . The decay indices 
of X-ray and optical afterglows are given with ax = {a - 3)/2 + 
(a - 1 l)(p - 2)/8 and a = {la - 25)/8 + (a - \l)(p - 2)/8 in the 
case of v m < vo < v c < vx, respectively. For the plateau phase of 
GRB 061 121, the X-ray and optical decay indices provide a con- 
sistent a within errors for possible values of p (2.2 > p > 2.8). 
In the case of p = 2.78, for example, we obtained a = 4.94± 
0.48 from ffxi and a = 5.05± 0.10 from aoi- For the shallow 
decay phase, assuming p = 2.68, for example, we calculated 
a = 5.99+ 0.44 from a X 2 and a = 5.23± 0.10 from a 02 . Thus, 



the decay indices of the plateau and shallow decay ph ases can be 
reprod uced with the prior activity model proposed by lloka et al.l 
(2006), while a temporal variation of a would be needed. 

Bo th models in iPanaitescu e t al.l (120061) and lloka et aP 
(l2006h can explain the observed light curves of GRB 061121, 
only when time variations in b, i, or a are allowed. The situation 
is further confusing when we consider the presence of the op- 
tical hump between the plateau and shallow decay phases. The 
discontinuity around the hump indicates a variation of the op- 
tical decay index, which means a further variation in b, i, or a 
during the hump. 

The hump in the optical light curve between the plateau and 
shallow decay phases is apparently not seen in the X-ray light 
curve. This implies that the dominant emitting regions are dif- 
ferent in X-ray and optical afterglows around the hump. The 
hump structure reminds us of the two-component jet model; the 
hump may be explained with a scenario that the emission from 
a narrow jet may dominate in the plateau pha se, while that from 
a wide jet beca me dominant after the hump dSheth et al] [2003: 
iPeng et a"Pl2005l) . In this case, the hump structure may appear 
when v m of the synchrotron emission from the wide jet passes 
the optical band. The classical synchrotron-shock model failed 
to reproduce the light curves possibly because of the composi- 
tion of the two components in the plateau and shallow decay 
phases. 

In the standard synchrotron-shock model, an increase of the 
density of the shock region would produce a hump in the opti- 
cal light curve which is n ot seen in X-rays if v c lies be tween the 
optical and X-ray bands ( Panaite scu and K umar 2000). Hence, 
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the optical hump of GRB 061 121 may also be reproduced if the 
external shock passed through a high density region in the inter- 
stellar medium at t = 2-6 ks. 

We finally note that there are several sources which ex- 
hibit optical li ght curves analogo us to GRB 0611 21, that is, 
GRB 021004 dUemura et al.| [2002b. GRB 50525A (iKlotz et al.1 
12001 lBlustin etal.1120061). G RB 060117 dJelmek et al.1 120061). 
GRB 060526 (iDai et al.ll2007l) . and GRB 061007 dMundell etal.1 
l2007h . For all of them, an early decay phase was terminated by a 
flattening or a hump at 10 2 ~ 4 s, which was followed by a steeper 
decay phase. An important point is that the decay indices before 
and after the flattening phase are different each other. This char- 
acteristic feature may commonly be observed in a group of GRB 
afterglows. If the emission during the early decay phase has a 
different nature from that during the later decay phase, the rela- 
tionship of those two components possibly causes the diversity 
in light curves of optical afterglows and the correlation between 
early X-ray and optical light curves. 

5. Summary 

We performed time-series photometry of the optical afterglow of 
GRB 061 121 with the 1.5-m Kanata telescope, and reported on 
a detailed study of the afterglow with published X-ray and opti- 
cal data. The decay index of the optical light curve was signif- 
icantly different between the plateau and shallow decay phases. 
The optical light curve possibly has a hump structure between 
the plateau and shallow decay phases, while no sign of such a 
hump is seen in the X-ray light curve. The different behavior in 
the optical and X-ray light curves indicates that they have dis- 
tinct emitting sources. The hump structure in the optical light 
curve may imply a passage of the typical frequency of the syn- 
chrotron emission from another forward shock distinct from the 
early afterglow. The optical decay index became same as the X- 
ray one in the late phase after the final break at ~ 4.6 x 10 4 s. 
In conjunction with the temporal evolution of SEDs, we propose 
that this break is caused by the passage of the cooling frequency 
at the optical band. In both the plateau and shallow decay phases, 
the observed decay and spectral indices are inconsistent with the 
standard synchrotron-shock model. They requires the variation 
of microphysical parameters in the shock region or the prior ac- 
tivity of the central engine. It is also possible that they are due to 
the composition of two forward shock components if the hump 
structure in the light curve was caused by another forward shock. 
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